Nonlinear transmission and light localization in photonic crystal waveguides 



Sergei F. Mingaleev and Yuri S. Kivshar 

Nonlinear Physics Group, Research School of Physical Sciences and Engineering 
Australian National University, Canberra ACT 0200, Australia 



We study the light transmission in two-dimensional photonic crystal waveguides with embedded 
nonlinear defects. First, we derive the effective discrete equations with long-range interaction for 
describing the waveguide modes, and demonstrate that they provide a highly accurate generaliza- 
tion of the familiar tight-binding models which are employed, e.g., for the study of the coupled- 
resonator optical waveguides. Using these equations, we investigate the properties of straight 
waveguides and waveguide bends with embedded nonlinear defects and demonstrate the possibility 
of the nonlinearity-induced bistable transmission. Additionally, we study localized modes in the 
waveguide bends and (linear and nonlinear) transmission of the bent waveguides and emphasize 
the role of evanescent modes in these phenomena. 



1. Introduction 



Photonic crystals are usually viewed as an optical ana- 
log of semiconductors that modify the properties of light 
similar to a microscopic atomic lattice that creates a 
semiconductor band-gap for electronall. One of the most 
promising applications of photonic crystals is a possibil- 
ity to create compact integrated optical deviceaa, which 
would be analogous to the integrated circuits in elec- 
tronics, but operating entirely with light. Replacing rel- 
atively slow electrons with photons as the carriers of in- 
formation, the speed and bandwidth of advanced com- 
munication systems can be dramatically increased, thus 
revolutionizing the telecommunication industry. 

To employ the high-tech potential of photonic crys- 
tals for optical applications and all-optical switching and 
waveguiding technologies, it is crucially important to 
achieve a dynamical tunability of their properties. For 
this purpose, several approaches have been suggested 
(see, e.g., Ref. ^). One of the most promising concepts 
is based on the idea to employ the properties of nonlin- 
ear photonic crystals, i.e. photonic crystals made from 
dielectric materials whose refractive index depends on 
the light intensity. Exploration of nonlinear properties 
of photonic band-gap materials is an important direc- 
tion of current research that opens a broad range of 
novel applications of photonic crystals for all-optical sig- 
nal processing and switching, allowing an effective way 
to create highly tunable band-gap structures operating 
entirely with light. 

One of the important concepts in the physics of pho- 
tonic crystals is related to the field localization on de- 
fects. In the solid-state physics, the idea of localization 
is associated with disorder that breaks the translational 



invariance of a crystal lattice and supports spatially lo- 
calized modes with the frequency outside the phonon 
bands. The similar concept is well-known in the physics 
of photonic crystals where an isolated defect (a region 
with different refractive index which breaks periodicity) 
is known to support a localized defect mode. An ar- 
ray of such defects creates a waveguide that allows a 
directed light transmission for the frequencies inside the 
band gap. Since the frequencies of the defect modes 
created by nonlinear dejects depend on the electric field 
intensity, such modes can be useful to control the light 
transmission. From the viewpoint of possible practical 
applications, spatially localized states in optics can be 
associated with different types of all-optical switching 
devices where light manipulates and controls light itself 
due to the varying input intensity. 

Nonlinear photonic crystals and photonic crystals with 
embedded nonlinear defects create an ideal environment 
for the observation of many of the nonlinear effects ear- 
lier predicted and studied in other branches of physics. 
In particular, the existence of nonlinear localized modes 
with the frequencies in the photonic band gaps has al- 
ready been predicted and demonstrated numerically for 
several models, of photonic crystals with the Kerr-type 
nonlinearityuQ. 

In this paper, we study the resonant light transmission 
and localization in the photonic crystal waveguides and 
bends with embedded nonlinear defects. For simplicity, 
we consider the case of a photonic crystal created by 
a square lattice of infinite dielectric rods, with waveg- 
uides made by removing some of the rods. Nonlinear 
properties of the waveguides are controlled by embed- 
ding the nonlinear defect rods. We demonstrate that 
the effective interaction in such waveguiding structures 
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is nonlocal, and we suggest a novel theoretical approach, 
based on the effective discrete equations, for describing 
both linear and nonlinear properties of such photonic- 
crystal waveguides and circuits, including the localized 
states at the waveguide bends. Additionally, we study 
the transmission of waveguide bends and emphasize the 
role of evanescent modes for the correct analysis of their 
properties. 

The paper is organized as follows. In Sec. 2, we intro- 
duce our model of a two-dimensional (2D) photonic crys- 
tal and provide a brief derivation of the effective discrete 
equations for the photonic-crystal waveguides created by 
removed or embedded rods, based on the Green function 
technique. In Sec. 3, we apply these discrete equations 
for the analysis of the transmission of straight waveg- 
uides, including the resonant transmission through an 
array of nonlinear defect rods embedded into a straight 
waveguide. Additionally, we discuss a link between our 
approach and the results obtained in the framework of 
the familiar tight-binding approximation often used in 
the solid-state physics models, and emphasize the im- 
portant role of the evanescent modes which can not be 
accounted for in the framework of the tight-binding ap- 
proximation. Section 4 is devoted to the study of waveg- 
uide bends. First, we discuss the localized modes sup- 
ported by a waveguide bend, and then we analyze the 
transmission of the waveguide bends. We demonstrate 
that the waveguide bends with embedded nonlinear de- 
fects can be used for a very effective control of light trans- 
mission. Section 5 concludes the paper with a summary 
of the results and further applications of our approach. 



2. Effective discrete equations 



radius = 0.18a and the dielectric constant Eq — 11.56 
(this corresponds to GaAs or Si at the wavelength ~ 1.55 



fj,m). For the electric field E{x,t) 



■ E{x\oj) po- 



larized parallel to the rods. Maxwell's equations reduce 
to the eigenvalue problem 



e{x) 



E{x\uj) = , 



(1) 



which can be solved by the plane-wave method^. A 
perfect photonic crystal of this type possesses a large 
(38%) complete band gap between uj = 0.303 x 27rc/a 
and to = 0.444 x 27rc/a (see Fig. |l|), and it has been ex- 
tensively employed during last few years for the study 
of bound stategj, transmission jiif light through shaqp 
bendsB'tj, waveguide-, branchesEJ and intersectionslld, 
channel drop filteisc^, nonlinear localized modes in 
straight waveguidesO and jdiscrete spatial solitons in per- 
fect 2D photonic crystals^. Recently, this type of pho- 
tonic crystal with a 90° bent waveguide has been fabri- 
cated in macro-porous silicon jadth a — 0.57 fim and a 
complete band gap at 1.55 /zttiE^. 
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In this Section, we suggest and describe a novel theoret- 
ical approach, based on the effective discrete equations, 
for describing many of the properties of the photonic- 
crystal waveguides and circuits, including the transmis- 
sion spectra of sharp waveguide bends. This is an im- 
portant part of our analysis because the properties of 
the photonic crystals and photonic crystal waveguides 
are usually studied by solving Maxwell's equations nu- 
merically, and such calculations are, generally speaking, 
time consuming. Moreover, the numerical solutions do 
not always provide a good physical insight. The effective 
discrete equations we derive below and employ further 
in the paper are somewhat analogous to the Kirchhoff 
equations for the electric circuits. However, in contrast 
to electronics, in photonic crystals both diffraction and 
interference become important, and thus the resulting 
equations involve the long-range interaction effects. 

We introduce our approach for a simple model of 2D 
photonic crystals consisting of infinitely long dielectric 
rods arranged in the form of a square lattice with the 
lattice spacing a. We study the light propagation in the 
plane normal to the rods, assuming that the rods have a 



Fig. 1. The band-gap structure of the photonic crystal 
created by a square lattice of dielectric rods with ro = 0.18a 
and eo = 11.56; the band gaps are hatched. The top right 
inset shows a cross-sectional view of the 2D photonic crystal. 
The bottom right inset shows the corresponding Brillouin 
zone, with the irreducible zone shaded. 

To create a waveguide circuit, we introduce a system 
of defects and assume, for simplicity, that the defects 
are identical rods of the radius located at the points 
a?„, where n is the index number of the defect rods. Im- 
portantly, the similar approach can be employed equally 
well for the study of the defects created by removing iso- 
lated rods in a perfect 2D lattice, and we demonstrate 
such examples below. 

In the photonic crystal with defects, the dielectric con- 
stant £(2^) can be presented as a sum of the periodic and 



defect-induced terms, i.e. £(2^) = 
Se{x) = ^ ed[E{x \ lj)] f{x 



:i{x) + Se{x), with 



(2) 



where f{x) — 1 for |x| < r^, and it vanishes otherwise. 



Equation (|l|) can be therefore written in the foUowing 
integral form 



(x|c.) - (^^)' / cPy G{x,y\uj)5e{y)E{y\u;) , 



(3) 



where G(a;, j/jo;) is the Green function of a perfect 2D 
photonic crystal (see, e.g., Ref. ||). 

A single defect rod is described by the function 
5e{x) — edf{x), and it can support one or more localized 
modes. Such localized modes are the eigenmodes of the 
discrete spectrum of the following eigenvalue problem. 



(x)=( 



c J 



d^y G{x,y\uji)sd.fiy)£i{y) , 



(4) 



where uji is the frequency (a discrete eigenvalue) of the 
Z-th eigenmode and £i{x) is the corresponding electric 
field. 

When we increase the number of defect rods (for 
example, in|-.cpder to create photonic-crystal waveg- 
uide circuitajiij) , the numerical solution of the inte- 
gral equation (^ becomes complicated and, moreover, 
it is severely restricted by the current computer facili- 
ties. Therefore, one of our major goals in this paper is 
to develop a new approximate physical model that would 
allow the application of fast numerical techniques, com- 
bined with a reasonably good accuracy, for the study of 
more complicated (linear and nonlinear) waveguide cir- 
cuits in photonic crystals. 

To achieve our goal, we consider the localized states 
created by a (in general, complex) system of defects (^) 
as a linear combination of the localized modes £i (x) sup- 
ported by isolated defects: 



Eix\u;)^J2'^nH^)£li^-Sn) 



(5) 



Substituting Eq. (||) into Eq. (||), multiplying it by 
£ii (x — Xn' ) and integrating with a?, we obtain a sys- 
tem of discrete equations for the amplitudes '0n ' of the 
l-th. eigenmodes localized at n-th defect rods: 



(6) 



where 



A,n = J (fx £i{x - Xn) £l'{x - Xn>) , 
f^in,mi^) = (■^) / £l'{x-Xn') X 



(fy G{x,y\uj) f{y-x„i)£i{y~Xn) . (7) 



It should be emphasized that the discrete equations (^- 
(0) are derived by using only the approximation provided 
by the ansatz (||). As can be demonstrated by compar- 
ing the approximate results with the direct numerical 
solutions of the Maxwell equations, this approximation 
is usually very accurate, and it can be used in many 
physical problems. 

However, the effective discrete equations (^ are still 
too complicated and, in some cases, they can be simpli- 
fied further still remaining very accurate. A good exam- 
ple is the case of the photonic crystal waveguides created 
by a sequence of largely separated defect rods. Such 
waveguides are known .as, the coupled-resonator optical 
waveguides (CROWs)Ejtj or coupled- cavity waveguides 
(CCWs)llII. For those cases, the localized modes are lo- 
cated at each of the defect sites being only weakly cou- 
pled with the similar neighboring modes. As is known, 
such a situation can be described very accurately by the 
so-called tight-binding approximation (see also Ref. 18). 



For our formalism, this means that ^" = /LtJ ^"^^ = 
for |n' — n| > 1 and \n' — m\ > 1. The most important 
feature of the CROW circuits is that thpE-bends are re- 
flectionless throughout the entire bandEJO. This is in 
a sharp contrast with the conventional photonic crys- 
tal waveguides created by a sequence of the removed or 
introduced defect rods (see e.g., Ref. || and references 
therein) where the 100% transmission through a waveg- 
uide bend is known to occur only at certain resonant fre- 
quencies. In spite of this visible advantage, the CROW 
structures have a very narrow guiding band and, as a re- 
sult, effectively they demonstrate a complete transmis- 
sion through the waveguide bend in a narrow frequency 
interval too. 

Below, we consider different types of photonic crystal 
waveguides and show that a very accurate simplifica- 
tion of Eqs. (^) is provided by accounting for an indi- 
rect coupling between the remote defect modes, caused 
by the slowly decaying Green function, fj,\ ^ for 
\n' — n\ < L, where the number L of effectively cou- 
pled defects lies usually between five and ten. As we 
show below, this type of interaction, which is neglected 
in the tight-binding approximation, is important for un- 
derstanding the transmission properties of the photonic- 
crystal waveguides. At the same time, we neglect a direct 
overlap between the nearest-neighhcui£igenmodes, which 
is often considered to be importantll3ilZl) , i.e. we consider 



A 



l.n 



5n,n' (with Si^ii being the Dirac delta func- 



tion) and ^\ — Q ioi n ^ m. Taking into account 
this interaction leads to negligible corrections only. 

Assuming that the defects support only the monopole 
eigenmodes (marked by I = 1), the coefficients (^ can be 
calculated reasonably accurately in approximation that 
the electric field remains constant inside the defect rods, 
i.e. £i{x) ~ f{x). This corresponds to the averaging 
of the electric field in the integtaLequation (|^) over the 
cross-section of the defect rodsE3El. In this case the re- 



suiting approximate discrete equation for the amplitudes 
of the electric fields En{ijj) = i'ni^) of the eigenmodes 
excited at the defect sites has the following matrix form: 



^Mn,m{L0)E,n{uj) = 
m 



(8) 



where J„.m(ijj) — ^{^^^{uj) is a coupling constant cal- 
culated in the approximation that £i{x) ^ f{x), so that 



J n. 



(9) 



is completely determined by the Green function of a per- 
fect 2D photonic crystal (see details in Rcfs. HJ^). 

First of all, we check the accuracy of our approxi- 
mate model (||) for the case of a single defect located 
at the point xq. In this case, Eq. (|^) yields the sim- 
ple result Jofl{ujd) — l/^d, and this expression defines 
the frequency Ud of the defect mode. In particular, ap- 
plying this result to the case when the defect is cre- 
ated by a single removed rod, we obtain the frequency 
LOd = 0.391 X l-Kcja which differs only by 1% from the 
value LOd — 0.387 x 2Tic/a calculated with the help of the 
MIT Photonic-Bands numerical codeQ. 



In Fig. H we plot the dispersion relation /ci(a;) calcu- 
lated by three different methods: first (solid curve) is 
calculated directly by the super-cell methodlD, and other 
two are found from Eq. (|l^) in the nearest-neighbor ap- 
proximation (L=l, dotted curve) corresponding to the 
tight-binding models, and also by taking into account 
the long-range interaction through the coupling between 
several neighbors {L=7, dashed curve). Indeed, we ob- 
serve a very good agreement for the conditions when the 
long-range interaction is taken into account. 
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Fig. 2. Dispersion relation for a 2D photonic-crystal 
waveguide (shown in the inset) calculated by the su- 
per-cell methodCI (solid), and from the approximate equations 
(0)-(|ll|) for L=7 (dashed) and L=l (dotted). The hatched 
areas are the projected band structure of a perfect 2D crystal. 



3. Straight waveguides 



A. Waveguide dispersion 



A simple single-mode waveguide can be created by re- 
moving a row of rods (see the inset in Fig. ||) . Assuming 
that the waveguide is straight (M„^m = Mn-m) and ne- 
glecting the coupling between the remote defect rods (i.e. 
Mn-m — for all \n — m\ > L), we rewrite Eq. (^) in 
the transfer-matrix form, Fn+i — TFn, where we intro- 
duce the vector Fn = { En, ■■• , En-2L+i } and the 
transfer matrix T — {Ti j} with the non-zero elements 



0.5 



Tij(w) 



Ml-, (a;) 



r„+i = l for J = 1,2, 



for i = l,2,. 
2L - 1 . 



.2L 



Solving the eigenvalue problem 
fiuj)$P = exp{ifcp(cj)}$P , 



(10) 



(11) 



we can find the 2L eigenmodes of the photonic-crystal 
waveguide. The eigenmodes with real wavenumbcrs 
kp{oj) correspond to the modes propagating along the 
waveguide. In the waveguide shown in Fig. H, there ex- 
ist only two such modes (we denote them as $^ and <I'^), 
propagating in the opposite directions {ki — — A;2 > 0). 
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Fig. 3. The same as in Fig. g but for two other types of 
waveguides better described by the tight-binding models. 



It is important to compare the results provided by our 
method with those obtained in the tight-binding approx- 
imation. For the waveguide shown in Fig. ^, such an 
approximation, rigorously speaking, is not valid, but its 
analog can be considered at least formally for the case of 
the nearest-neighbor interaction when in Eqs. ([lC|)-(|ll|) 
we take L=l. The interaction between the remote rods 
cannot be neglected as soon as we study the waveguides 
created by removing (or inserting) all rods along the row 
or more complicated structures of this type. In such a 
case, as is seen from Fig. |[ the dispersion relation found 
in the tight-binding approximation is incorrect and, in 
order to obtain accurate results, one should take into 
account the coupling between several defect rods. We 
verify that this statement is also valid for multi-mode 
waveguides, e.g. those created by removing several rows 
of rods. 
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Fig. 4. Transmission coefficients of an array of nonlinear 
defect rods calculated from Eqs. (^)-(^) with L = 7 in the 
linear limit of very small |atp (solid curves) and for the non- 
linear transmission when the output intensity is |Qtp = 0.25 
(dashed curves), for different numbers of the defects. We use 



nonlinear defect rods with the dielectric constant e 
they are marked by open circles on the right insets. 



(0) 



7; 



However, for the waveguides of a different geometry, 
when only the second (or third, etc) rods are removed, 
all methods provide a reasonably good agreement with 
the direct numerical results, as is shown for two exam- 
ples in Fig. 1^. In this case, the waveguides are cre- 
ated by an array of cavity modes, and they are simi- 



lar to the. CROW structures earlier analyzed by several 
authorsliErEZl. Thus, the dispersion properties of CROWs 
(or similar waveguides) can be described with a good ac- 
curacy by the tight-binding approximation; our new ap- 
proach confirms this conclusion, and it provides a simple 
method for the derivation of the approximate equations. 



B. Resonant transmission of an array of defects 

In addition to the propagating guided modes, in photonic 
crystal waveguides there always exist evanescent modes 
with imaginary kp. These modes cannot be accounted for 
in the framework of the tight-binding theory that relies 
on the nearest-neighbor interaction between the defect 
rods. Although the evanescent modes remain somewhat 
"hidden" in straight waveguides, they become crucially 
important in more elaborated structures such as waveg- 
uides with embedded linear or nonlinear defects, waveg- 
uide bends and branches. In these cases the evanescent 
modes manifest themselves in several different ways. In 
particular, they determine non-trivial transmission prop- 
erties of the photonic-crystal circuits considered below. 

As the first example of the application of our approach, 
we study the transmission of a straight waveguide with 
embedded nonlinear defects. Such a structure can be 
considered as two semi-infinite straight waveguides cou- 
pled by a finite region of defects located between them. 
The coupling region may include both linear (as a do- 
main of a perfect waveguide) and nonlinear (embedded) 
defects. We assume that the defect rods inside the cou- 
pling region are characterized by the index that runs 
from a to &, and the amplitudes Em (m = a, . . . ,b) of 
the electric field at the defects are all unknown. We num- 
ber the guided modes (|l^) in the following way: p = 1 
corresponds to the mode propagating in the direction of 
the nonlinear section (for both ends of the waveguide), 
p — 2 corresponds to the mode, propagating in the oppo- 
site direction, p = 3, L -I- 1 correspond to the evanes- 
cent modes which grow in the direction of the nonlinear 
section, and p = L -\- 2, ...,2L correspond to the evanes- 
cent modes which decay in the direction of the nonlinear 
section. Then, we can write the incoming and outgoing 
waves in the semi-infinite waveguide sections as a super- 
position of the guided modes: 



L+1 

E 

p=3 



(12) 



for m = a — 2L, a — 1, and 

L + l 



E°f = at$2 



rn—o ' / J r'p 
p=3 



P 

m — b ' 



(13) 



for TO = 6+1, b+2L, where /3™ and are unknown 
amplitudes of the evanescent modes growing in the di- 
rection of the nonlinear section, whereas a,, at and Ur 



are unknown amplitudes of the incoming, transmitted, 
and reflected propagating waves. We take into account 
that the evanescent modes with p > L + 1 (growing in 
the direction of waveguide ends) must vanish. Now, sub- 
stituting Eqs. into Eq. (|), we obtain a system 
of linear (or nonlinear, for nonlinear defects) equations 
with 2L + b — a + l unknown. Solving this system, we find 
the transmission coefficient, T = \at/ai\'^, and reflection 
coefficient, R = jar/aip, as functions of the light fre- 
quency, Lu, andintensity, la^p or |atp. Recently, we have 
demonstratedcO that the linear transmission properties 
of the waveguide bends are described very accurately by 
this approach. Below, we study nonlinear transmission 
of the photonic-crystal waveguides and waveguide bends. 
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Fig. 5. Bistability in the nonlinear transmission of an 
array of five nonlinear defect rods shown in Fig. ^(b). 

In Fig. ^ we present our results for the transmission 
spectra of the straight waveguides (created by a row of 
removed rods) with an array of embedded nonlinear de- 
fects. We assume throughout the paper that all non- 
linear defect rods are identical, with the radius — 

ro = 0.18a and the dielectric constant, Sd — e?^ -I- |i?nP 



(with 



-(0) 



7), which grows linearly with the light in- 



tensity (the so-called Kerr effect). In the linear limit, 
the embedded defects behave like an effective resonant 
filter, and only the waves with some specific resonance 
frequencies can be effectively transmitted through the 
defect section. The resonances appear due to the exci- 
tation of cavity modes inside the defect region, whereas 
a single defect does not demonstrate any resonant be- 
havior. When the intensity of the input wave grows, the 
resonant frequencies found in the linear limit get shifted 
to lower values. The sensitivity of different resonances 
to the change of the light intensity is quite different and 
may be tuned by matching of the defect parameters. The 
nonlinear resonant transmission is found to possess bista- 
bility, similar to another problem of the nonlinear trans- 
mission (see, e.g., Refs. ^l|). The bistable transmission 
occurs for the frequencies smaller then the resonant, in 
a linear limit, frequency (see Fig. |^). 



C. An optical diode 



An all-optical "diode" is a spatially nonreciprocal device 
which allows unidirectional propagation of a signal at 
a given wavelength. In the ideal case, the diode trans- 
mission is 100% in the "forward" propagation, while it 
is much smaller or vanishes for "backward" (opposite) 
propagation, yielding a unitary contrast. 

The first study of the operational mechanism for a pas- 
sive optical diode based on a photcmip_band gap mate- 
rial was carried out by Scalora et a/.E3o. These authors 
considered the pulse propagation near the band edge of a 
one-dimensional photonic crystal structure with a spatial 
gradiation in the linear refractive index, together with a 
nonlinear medium response, and found that such a struc- 
ture can result in unidirectional pulse propagation. 
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Fig. 6. Transmission coefficients of an asymmetric array 
of nonlinear defect rods calculated for the same parameters 
as in Fig. ^. 



0.15 




Incoming light, I a. I 



Fig. 7. Nonlinear transmission of the optical diode for 
the forward (see top of Fig. ^ and backward (see bottom of 
Fig. ^) directions at the light frequency lo — 0.326 x 2-Kc/a. 

To implement this concept for the waveguide geometry 
discussed above, we consider the asymmetric structure 
made of four nonlinear defect rods, as shown in the right 
insets on Fig. O. Figure shows the transmission spectra 



of such an asymmetric structure in the opposite direc- 
tions indicated by two arrows in the right insets. As is 
seen, in the hnear limit the transmission is characterized 
by two resonant frequencies and does not depend on the 
propagation direction. However, since the sensitivity of 
both resonant frequencies to the change of the hght in- 
tensity is different for the "forward" (see Fig. ^ top) and 
"backward" (see Fig. ||, bottom) propagation directions, 
the transmission becomes, in the vicinity of resonant fre- 
quencies, highly asymmetric for large input intensities. 
This results into nearly unidirectional waveguide trans- 
mission, as shown in Fig. ^ 

In contrast to the perfect resonators used in Fig. ^, the 
transmission of the asymmetric structure under consid- 
eration is not very efficient at the resonant frequencies. 
However, we believe that the optical diode effect, with 
much better efficiency, can be found in other types of 
the waveguide geometry and a unitary contrast can be 
achieved by a proper optimization of the waveguide and 
defect parameters, that can be carried out by employing 
our method and the effective discrete equations derived 
above. 



4. Waveguide bends 



In addition to the non-trivial transmission properties, 
the evanescent modes manifest themselves in creation of 
localized bound states in the vicinity of the waveguide 
bends or branches. To demonstrate this effect, we con- 
sider the simplest case of a sharp bend of the waveguide 
created by a row of removed rods. As was shown in 
Ref. ^, in the cases when the waveguide bend can be 
considered as a finite section of a waveguide of differ- 
ent type, the bound states correspond closely to cavity 
modes excited in this finite section. However, such a sim- 
plified one-dimensional model does not describe correctly 
more complicated casesQ, even the properties of the bent 
waveguide depicted in Fig. ^. The situation becomes 
even more complicated for the waveguide branchestj. In 
contrast, solving the discrete equations (^ we can find 
the frequencies and profiles of the bound states excited 
in an arbitrary complex set of defects. For particular 
case of the waveguide bend shown in Fig. ||, we find two 
bound states localized at the bend and their profiles (cf. 
our Fig. H with Fig. 9 in Ref. ||). It should be noted that 
the frequencies of the modes are found from Eq. (H) with 
the accuracy of 1.5%. 

Additionally, the evanescent modes determine the 
non-trivial transmission properties of the waveguide 
bends which can also be calculated with the help of 
our discrete equations. To demonstrate these features, 
we consider a bent waveguide consisting of two coupled 
semi-infinite straight waveguides with a finite section of 
defects between them. The finite section includes a bend 
with a safety margin of the straight waveguide at both 
ends. 
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Fig. 8. Electric field En for two bound states supported 
by a 90° waveguide bend shown in the top. Center of the 
bend is located at n = 0. 
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Fig. 9. Reflection coefiicients calculated by the fi- 
nite-difference time-domain method (dashed, from Ref. ^) 
and from Eqs. (||)-(|l|) with L = 7 (full lines) and L = 2 
(dotted, only in the top plot), for different bend geometries. 



Similar to what we discussed in Sec. 3B for the 
straight waveguides, we solve the system of effective dis- 
crete equations to find the transmission, T = \at/ai\'^, 
and reflection, R — |ar/aip, coefficients of the waveg- 
uide bends. In Fig. ^ we present our results for the 
transmission spectra of several types of bent waveguides. 



which have been discussed in Ref. ^, where the possibihty 
of high transmission through sharp bends in photonic- 
crystal waveguides was first demonstrated. Wc com- 
pare the reflection coefficients calculated by the finite- 
difference time-domain method in Ref. 1 (dashed lines) 
with our results, calculated from Eqs. (|)-(13) for L — 7 
(full lines) and for L — 2 (dotted line in the top plot). 
As is clearly seen, Eqs. (|)-(|l|) provide a very accu- 
rate method for calculating the transmission spectra of 
the waveguide bends, if only we account for long-range 
interactions, ft should be emphasized that the approxi- 
mation, in which the only next-neighbor interactions are 
taken in to account, is usually too crude, while the tight- 
binding theory incorrectly predicts a perfect transmis- 
sion for all guiding frequencies. 
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Fig. 10. Transmission of a waveguide bend with three 
embedded nonlinear defect rods in the linear (solid curve) 
and nonhnear (dashed curve) regimes. Defect rods have the 
dielectric constant e^"' = 7, and they are marked by open 
circles. 
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Fig. 11. Bistable nonlinear transmission through the 
waveguide bend shown in Fig. |l^, for the light frequency 
Lu = 0.351 X 2Tvc/a. 



The resonant transmission can be modified dramat- 
ically by introducing both linear and/or nonlinear de- 
fects into the waveguide bends. To illustrate such fea- 
tures, we consider the waveguide bend with three em- 
bedded nonlinear defects, as is depicted in Fig. |l^ on 
the right inset, where these defects are shown by open 
circles. In the linear regime, such a sharp bend behaves 
as an optical threshold device that efficiently transmits 
the guided waves with frequencies above the threshold 
one, but completely reflects the waves with the lower 
frequencies. The transmission coefficient of this waveg- 
uide bend in the linear limit is shown in Fig. by a solid 



curve. When the input intensity increases, the threshold 
frequency decreases, extending the transmission region 
(see a dashed curve in Fig. |l^). The resulting transmis- 
sion as a function of the input intensity demonstrates a 
sharp nonlinear threshold character with an extremely 
low transmission of the waves below a certain (rather 
small) threshold intensity, see Fig. |ll|. 



5. Conclusions 



We have suggested a novel theoretical approach for de- 
scribing a broad range of transmission properties of the 
photonic crystal waveguides and circuits. Our approach 
is based on the analysis of the effective discrete equations 
derived with the help of the Green function technique, 
and it generalizes the familiar tight-binding approxima- 
tions usually employed to study the coupled-resonator or 
coupled-cavity optical waveguides. The effective discrete 
equations we have introduced in this paper emphasize 
the important role played by the evanescent modes in 
the transmission characteristics of the photonic crystal 
circuits with waveguide bends and/or embedded defects. 
Employing this technique, we have studied the proper- 
ties of several important elements of the (linear and non- 
linear) photonic crystal circuits, including a nonlinear 
bistable transmitter and an optical diode created by an 
asymmetric structure of nonlinear defects. We believe 
that our approach can be useful for solving more com- 
plicated problems, and it can be applied to study the 
transmission characteristics of the waveguide branches, 
channel drop filters, etc. 

The work has been partially supported by the Aus- 
tralian Research Council. 
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